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Abstract 

We investigate the Yang-Lee edge singularity on non-planar random graphs, which we consider 
as the Feynman Diagrams of various d = field theories, in order to determine the value of the edge 
exponent a. 

We consider the hard dimer model on 0^ and random graphs to test the universality of the 
exponent with respect to coordination number, and the Ising model in an external field to test its 
temperature independence. The results here for generic ( "thin" ) random graphs provide an interesting 
counterpoint to the discussion by Staudacher of these models on planar random graphs 



1 Introduction 



The work of Yang and Lee [|[ ^ , later expanded by various other authors , on the behaviour of spin 
models in complex external fields has provided an important paradigm for the understanding of the 
nature of phase transitions. In brief, Yang and Lee observed that the partition function of a system 
above its critical temperature Tc was non-zero throughout some neighbourhood of the real axis in the 
complex external field plane. As T — > Tc+ the endpoints of loci of zeroes moved in to pinch the real axis, 
signalling the transition. When such endpoints occur at non-physical (i.e. complex) external field values 
they can be considered as ordinary critical points with an associated edge critical exponent. The picture 
was later extended by Fisher to temperature driven transitions by considering the analyticity properties 
of the free energy in the complex temperature plane 

A few equations to flesh this out would perhaps not go amiss. On a finite graph G„ with n vertices 
the free energy of an Ising-like spin model can be written as 

n 

F{Gn,f3,z) = -nh-lnY[{z-Zk{f3)) (1) 

fe=i 

where the fugacity z = cxp(— 2/i), and h is the (possibly complex) external field. The Zk{(3) are the 
Yang-Lee zeroes, which in the thermodynamic limit form dense sets on curves in the complex z plane. In 
the infinite volume limit rt — > oo the free energy per spin is 

F(Goo, /3, z) ^-h- f d0p{l3, 9) Iniz - e'^) (2) 



where p(/3, 9) is the density of the zeroes, which can be shown to appear on the unit circle in the complex 
z plane in the Ising case (the Yang-Lee circle theorem). For T > Tc or, if one prefers /3 < /3c, there is a 
gap with p{f3, 9) ~ for \9\ < 9o, and at these edge singularities we have 

p{f3,9)^{9~0or (3) 

which defines the Yang-Lee edge exponent cr. This also implies M ^ {9 — 9o)'^ . Various finite size scaling 
relations relate the Yang-Lee exponent to the other critical exponents and can be used in numerical 
determinations of critical behaviour M . 

The Yang-Lee circle theorem of gTH] guarantees that the roots of the partition function of the Ising 
model on a fixed graph Gn lie on the unit circle, but it does not guarantee that this should be the case 
when the partition function is defined by a sum over some class of random graphs for each n 

Z^ = Y,Z{Gn) (4) 

where Z(Gn) is the partition function on a given graph in the class. This is the case for models of 2D 
quantum gravity, where the sum is over planar (j)"^ , (j)'^ . . . random graphs, and in this paper where we will 
will consider a sum over thin, non-planar random graphs |^. The work of Staudacher showed that, 
nonetheless, the zeroes did appear on the unit circle for planar graphs, which has recently been confirmed 
by numerical finite-size scaling investigations using both series expansions and Monte-Carlo simulations 
by Ambj0rn et.al. 

In this paper we will consider a partition function of the form equ.(^ for thin random graphs. Spin 
models on such thin random graphs are of interest because they give a way of investigating mean-field 
efl[ects thanks to their tree-like local structure ^ . The advantage of using the thin random graphs in such 
investigations over genuine tree-like structures such as the Bethe lattice is that boundary effects are 
absent. The complications of being forced to consider only sites deep within the lattice which occur on 
the Bethe lattice are thus absent. The motivation for the current work is the calculation of the exponent 

^We shall call such graphs "thin" random graphs throughout this paper as they appear as the scalar limit of the matrix 
fatgraphs that are relevant for discussions of 2D gravity. 
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X iV„ = — -—^ I J^^,^ exp(-5), (5) 



a for thin random graphs and the testing of its universahty. The calculation of the edge exponent in the 
Ising model context also allows us to check its temperature independence. A mean-field value for a (i.e. 
1/2) would demonstrate that the mean-field nature of critical behaviour on such random graphs models 
extended to complex couplings. 

We shall use the approach of |l2[ |l^ to solve both the hard dimer model and the Ising model itself in 
a complex external field. The requisite ensemble of thin random graphs is generated by considering the 
scalar limit of a matrix model. In all cases the partition function of equ.(^) is given by an integral of the 
form 

27ri / A2"+i J 2TTVd^ 
where K is the inverse of the quadratic part of the action S, the (j)i are the fields which give the appropriate 
decoration of the graph, and A is the vertex coupling. The factor iV„ counts the number of undecorated 
graphs in the class of interest, and generically grows factorially with n. A given graph appears as a 
particular "Feynman diagram" in the expansion of equ.(^ and the integration over A picks out graphs 
with 2n vertices. The coupling A is irrelevant for the discussion of critical behaviour as it may be scaled 
out of the action and hence any saddle point equations. In the large n limit the integral in equ.(^ may 
be evaluated by saddle point methods. Phase transitions appear when an exchange of dominant saddle 
points occurs, either continuously giving a second order transition, or discontinuously giving a first order 
transition. The saddle point integrals which appear are the o? = equivalent of those in instanton and 
large orders calculations in field theory |l^ ^ . As we are taking an n ^ oo limit at the start of our 
calculations, we are unable to explicitly calculate zeroes on finite lattices and verify the Yang-Lee circle 
theorem. However, we shall see that the end-points of the loci of zeroes do lie on the unit circle in the 
complex fugacity plane. 

One important property of the Yang-Lee edge exponent a in all the models examined so far is that it 
is independent of f3 (for f3 < f3c). It is therefore possible in general to get a quick and dirty determination 
of its value for Ising models by taking the so-called hard dimer limit, which corresponds to (3 ^ 0, h ^ 
and has certain simplifying features compared with the general case. We will calculate the edge exponent 
for the hard dimer model on (p^ and 0"* random graphs in the next section. A calculation of the exponent 
for the Ising model proper in the section which follows this will investigate whether the temperature 
independence still holds. 

2 Hard Dimer Models 

The partition function of the hard dimer model on a given graph is defined [0 by 

e(G„) 

e(G„) = l+ ^ 0„(z)C (6) 

where C is a dimer activity and 0„(i) is the number of ways of placing i dimers on the e(G„) edges of 
the graph Gn such that at most one dimer is attached to each vertex (hence "hard"). Precisely this 
expression appears in taking the limit (3 0,h ^ in the high temperature series for the Ising model 
partition function. The role of h in the definition of the edge singularity exponent is taken by ^ and we 
have 

f~(c-c.r (T) 

where Q is the appropriate n oo limit of 0(G„). For both planar and thin random graphs we are 
interested in a sum 

e„ = ^e(G„) (8) 

and in Q the appropriate two matrix integral to generate the dimer partition function on 0'^ and 0'* 
planar graphs was written down. For thin 0^ graphs the required action for insertion into equ.(|^) is 

S=\{x' + y')-lx^-VCyx'. (9) 
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where the x propagators represent the unoccupied hnks and the y propagators the dimers. The two 
vertices are shown in Fig. 1 . The -y/C weight appears because each end of the dimer contributes a ^/Qyx^ 
vertex. We have scaled out A for clarity. Similarly for (f)^ graphs we find 

S=\{x' + y')-\x'-^Cyx\ (10) 

The saddle point equations dS/dx ~ dS/dy — for both equs.(p|,|lO|) can be easily solved. Concentrating 
on the (f)^ case for simplicity, we find 



. = -i^^™^ y=^—^ ill) 
so the resulting saddle point action is 

s^^^(i2C-sc.JTTk + '^K^ + i-VTTs^). (12) 

We can see that the saddle point solution presents a singularity at a negative value, Co = —1/8. The free 
energy is given by the logarithm of the action to leading order in l/n, so we would expect an inverse 
square root divergence at Co when we differentiate In S twice if a were equal to its mean-field value of 
1/2. This is, indeed, the case. Writing C, = — 1/8 + e and expanding we find 

d^e InS 96V2 

As y appears at most quadratically in equ. (^) an alternative approach Q is to integrate it out to get the 
action 

S^l{x'+y')-lx^-lCx' (14) 

which has the saddle point solution 



(15) 

and displays an identical divergence ~ 96\/2/-ye at Co = —1/8 to equ. (|l^) . The geometrical picture of 
such an integration on the y variables is that all the dimers are collapsed to give new vertices and 
assigned the appropriate weight. Whichever way the calculation is carried out, the appearance of the 
square root divergence confirms that ct = 1/2, which is the mean-field value of the edge exponent. 

The universality of the result with respect to the coordination number of the vertices can be confirmed 
by solving the saddle point equations for equ.(p^), which we do not reproduce here, or by integrating out 
y to give 

S-l{^' + y')-\x'-lcx' (16) 

which has the saddle point solution 

. ^ + (17) 

V6C 

leading to the saddle point action 



1 (-1 + V1 + 12C)(24C + 1- V1 + 12C) 
432 C^ 



s^_±_ y-^^^^^^-^jyy^^^^-^^^^-^\ (18) 



When expanded around the singularity at Co = —1/12 it also gives a square root divergence 

d^e InS 3673 

^Also taken in jlj for the planar case in a matrix model calculation. 



(19) 
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The exponent a is thus seen to be independent of the coordination number of the random graphs (three 
and four for (j)^ and 0^ graphs, respectively) on which the dimers are placed. 

The hard dinier calculation represents a determination of the edge exponent at one point on the (/i, (3) 
plane. We now turn to the solution of the Ising model in an external field in order to confirm this value 
for generic points on the singular h{[3) line. 



3 The Ising Model in an External Field 

The action for the Ising model on (jy^ graphs in an external field may be written as 

S^\{^'+ y') - cxy - le'^x^ - \e-S' (20) 

where c = exp(— 2/3). The transition point in the model is determined by first solving the saddle point 
equations dS/dx — dS/ dy — and then using these solutions to determine at which point the Hessian of 
the second partial derivatives is zero pT[ |. This will pick up any continuous transitions that are present. 
The net result of these (lengthy) calculations is the following formula for h(c), the curve in the h, c plane 
along which the Hessian is zero 



'2c (1 + 18c2 - 27c4 ± (1 - 9c2) Vl - lOc^ + Oc^ ) 
exp(Mc)) - ±^^- (21) 

It is perhaps worthwhile to consider the solution in zero external field at this point for orientational 
purposes. In that case we have at high temperatures 

a; y = 1 - c (22) 

which bifurcates in a mean-field magnetisation transition at c = 1/3 to the low temperature solutions 



1 + c + Vl - 2c - 3c2 



1 + c- Vl-2c-3c2 

y = 2 ■ ^ ' 

valid for c < 1/3. The distinguished role of the zero-field critical point c = 1/3 is clear in equ.(^l]). 
exp(ft,(c)) develops an imaginary part for c > 1/3 (in the high temperature phase) and if we plot its 
modulus as in Fig. 2 we can see clearly that | exp(/i(c))| — 1 for c > 1/3. This shows that the endpoints, 
at least, of the line of zeroes lie on the unit circle. Without explicit finite size calculations, or simulations 
of the model in a complex field in the style of ||^ we cannot say that the zeroes lie on the unit circle, but 
it is reasonable to conjecture that they do, just as for planar graphs. 

Extracting the critical exponent from the magnetisation turns out to be easiest proposition for the 
Ising model in a field. The expression for the magnetisation in this case is 

"<«(")) - s^Tl^ (21) 

and it is directly related to the density of zeroes on the unit circle by 

p{e) - lim Re M (re'^) . (25) 

We proceed by substituting the saddle point solutions for x,y in field that led to equ.(Elh into equ.(p4) 



before taking the limit in equ.(25) above. The positions of the critical endpoints on the unit circle can 



be extracted by examining the discontinuities in the resulting expression, and are given by 

. / N 1 1 / (9c^ - l)\/l-10c2 + 9c4\ 

6*0 c = ± - tan^i ^ 26) 

' 2 I 1 + 18c2 - 27c4 j ^ ' 
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which is also consistent with equ.(pTl) for exp(/i). The endpoints can be seen to move in to pinch the real 
axis, ^0 0± as c ^ 1/3+ , confirming the Yang-Lee picture of the transition. 

An expansion of the density of zeroes around 9q for generic c still rapidly degenerates into considerable, 
and not very illuminating, algebraic complexity. However, fixing c and examining various points along 
the critical line in equ.(pl|) reduces the level of difficulty to that of the hard dimer calculations in the 
previous section. For any given value of c > 1/3, we do indeed find that p{9) ^ {9 ~ 6*0 (c))^/^ by using 
equ.(25). In Fig. 3 the density of zeroes is plotted against 9 — 9o{c) for c = 1/2 and the square root nature 
of the singularity is evident, as can be confirmed by fitting to the curve for this and other values of c. 



4 Discussion 

Both the hard dimer calculations and those for the full Ising model in an external field produce the 
mean-field value for the edge exponent. Given the body of previous results in |pl| showing mean-field 
behaviour in various models on random graphs (and the observation in the first of |pl]| that the saddle 
point equations in the mean-field models were identical in content to the recursion relations use to solve 
the models on trees) this is no great surprise, though it does provide the first confirmation that the 
mean-field critical behaviour on random graphs extends to critical phenomena at complex couplings. 
The inherent simplicity of the saddle point equations for random graphs made obtaining expressions for 
the singular behaviour and the position of the critical endpoints a simple task in the hard dimer model 
and also allowed for a demonstration of universality by examining both (f)^ and (f>^ graphs. The solution 
of the Ising model in field was rather more forbidding, and we have not presented many of the resulting 
elephantine equations here, but it was still possible to give a simple demonstration by semi-numerical 
means that a = 1/2 along the critical curve. 

We have also seen that for the Ising model the singular endpoints of the line of zeroes stay on the 
unit circle in the complex fugacity plane, which provides support for the Yang-Lee circle theorem with 
partition functions of the form in equ.(^ that involve sums over thin random graphs. This naturally 
leads to the conjecture that all the Yang-Lee zeroes lie on the unit circle for the Ising model on thin 
random graphs, just as on planar graphs, which could be confirmed by a finite size scaling analysis for 
the thin random graph model in the manner of that carried out in |^ for planar graphs. Another possible 
extension of the present work would be to consider the Fisher zeroes in the model by examining the 
behaviour in the complex temperature plane. 

Indeed, a more comprehensive investigation of complex phases for various spin and vertex models on 
both planar and non-planar random graphs, in the manner of that carried out by Matveev and Shrock [p^ 
for regular lattices, might prove illuminating. In particular, it would be interesting to see if the complex 
temperature and complex field singularities with atypical (and lattice dependent) exponents found in [p^ 
had their counterparts in random graph models. 
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(a) (b) 



Figure 1: The two vertices in the cf)^ dimer model, with the dimer edge drawn in bold, (a) carries a 
weight of A\/C and (b) carries a weight of A/3. 
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Figure 2: The modulus of exp(ft,) along the Yang-Lee transition line for c > 1/3 can be seen to be one 
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Figure 3: The density of zeroes p (suitably scaled for plotting convenience) is plotted against e = 6 — d{c) 
for c = 1/2 to show the square root edge singularity. 
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